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Abstract — This paper presents the value-at-risk (VaR) analysis of J.C. Penney Company Inc. (JCP) stock daily negative log 
returns between 1993 and 2018. The statistical properties of JCP are thoroughly examined and a series of diagnostic tests 
are conducted to check the conditions of the time series data over the two decades. The GARCH and EGARCH models with 
normal distribution and Student’s t-distribution are used to estimate the volatility and VaR of the stock. By analyzing VaR, 
we show that there is currently a high risk of investing in JCP stock. In addition, this paper examines the calendar effects 
and seasonality of JCP stock through the fundamental properties of the data as well as the VaR. We compare the 
performance of the stock in four quarters which further confirms our result that JCP stock is at immense risk at this point in 
time. These results are valuable for anyone interested in evaluating and forecasting JCP stock. The methodology we use is 
applicable to any other stock that meets our test conditions and is more accurate and realistic in predicting volatility and 
VaR than the commonly used standard normal distribution based VaR model. 
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I. Introduction 

Over the past two decades, the proliferation of the Internet has contributed to the development of e-commerce and has 
affected the market share of traditional retailers. Purchasing goods online as a convenient alternative to in-store shopping 
largely saves people’s time and is increasingly popular among customers. J.C. Penney Company Inc. is a US department 
store chain that operates more than 1000 stores across the United States. It was founded by James Cash Penney and William 
Henry McManus [1] in 1902. Most J.C. Penney stores are located in shopping centers, and their business mainly includes 
sales of clothing, cosmetics, household items, jewelry and cookware. As J.C. Penney is one of the largest apparel and home 
retailers in the United States, investors are deeply concerned about the performance of J.C. Penney Company Inc. stock 
(NYSE: JCP). In this paper, we will investigate the risk of buying JCP stock as well as the calendar effects of the JCP stock 
returns. We want to provide investors with useful investment advice. 

In 2017, Caroline, Emma, Madelon, Mikkel and Marc from Columbia Business School conducted a research project on J.C. 
Penney named “ Competingfor Survival: A Turnaround of Department Store J.C. Penney” [2]. They analyzed the company’s 
performance from an operational perspective and provided a variety of strategies to help the company better organized. So 
far, we have found a lack of statistical analysis of J.C. Penney's paper. In this paper, we conducted a complete time series 
analysis of the JCP stock returns. We focus on the statistical characteristics of JCP stock and draw conclusions only from the 
statistics. 

Risk management is a key process for making investment decisions. In order to control risk in an investment, we need to first 
determine the amount of risk involved in the investment, and then decide to either accept or alleviate the risk. Standard 
deviation, beta, value at risk (VaR) and expected shortfall are common measures to quantify the risk. In this paper, we use 
VaR as the primary tool for measuring the risk of buying stocks. Value-at-risk is a statistical indicator of the riskiness of 
financial entities or portfolios. It is defined as the maximum dollar amount that is expected to be lost at a predetermined 
confidence level for a given time frame. The stock market crash in 1987 triggered the innovation of VaR. It was developed as 
a systematic approach to separating extreme events from daily price changes. In 1994, it was extended by J.P. Morgan who 
launched the Risk Metrics and published the methodology [3]. Nowadays, VaR has become one of the most commonly used 
measures of market risk in the financial industry. 

In order to calculate the VaR of the stock, we need to accurately predict the price of the stock. The main characteristic of a 
stock is its return. We model the negative log return series to derive estimates of volatility and VaR. It is well known that 
financial markets are highly volatile and the periods of high volatility tend to persist for some time before the market returns 
to a more stable environment (Tsay,2005)[4j. The autoregressive method helps to build a more accurate and reliable 
volatility model. The Autoregressive Conditional Heteroskedasticity (ARCH) model was originally introduced by Engle 
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(1982)[5]. ARCH model and its extensions such as GARCH (Bollerslov, 1986)[6] and EGARCH (Nelson, 1991) [7] are 
among the most popular models for forecasting market returns and volatility. See for example Vesna Bucevska (2013)[8], 
Zhen Yao Wong et al. (2016)[9], Julija Cerovic Smolovic et al. (2017)[10] and Mahsa Gorji, Rasoul Sajjad (2017) [11], who 
derived VaR estimation from GARCH - type models. 

In terms of the organization of the paper, we analyze the basic statistical properties of the time series of negative log returns 
for JCP in Section 2, and test for the normality and autocorrelations of the series. In Section 3, we fit the data with GARCH 
and EGARCH models and estimate the VaR of the data. In Section 4, we group the data by quarter and perform time series 
analysis within the group, and then we analyze the results for each group to find calendar effects 

II. Data Description and Statistical Test Results 


2.1 Data Description 

In this paper, we examine the daily JCP stock price time series for a 25-year period. There were 6431 data points from 
January 29, 1993 to August 10, 2018. All data comes from Yahoo Finance. Let P t denote the daily adjusted closing price of a 

stock, where t is an integer representing the day. We use the negative log return L f to characterize the stock price time 
series. L f is defined as: 


L ( = -100• log2L = I00(log P t ] - log P t ) 

Pf-l 


( 1 ) 
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Figure 1: The upper panel is adjusted closing price for JCP; the lower panel is daily negative log returns 
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JCP.rtn [1993-02-01/2018-08-10] 
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FIGURE 2: The daily squared log returns for JCP 


The lower panel in Fig. 1 shows the time series plot of daily negative log returns for JCP stock. We can see from the graph 
that roughly from 2013 to 2018, there are more peaks in the upper half of the panel. Since we model relative loss instead of 
relative return, more positive peaks indicate more positive relative loss of the stock. Fig. 2 shows the time series plot of daily 
squared log returns of JCP stock. It can be seen that since 2011, the stock has become more volatile than before, and the 
stock volume shown in Fig. 1 has become higher after 2011. There are a lot of documented evidence of high stock trading 
volume closely related to the volatility of returns; see for example Barron, Orie E., David G. Harris, and Mary Stanford 
(2005[12]), K. Ravichandran, Sanjoy Bose (2012 [13]), Andrey Kudryavtsev (2017[14]). 

From the price chart in Fig. 1, JCP stock reached its highest price in 2007 and then fell sharply in the following year. From 
2012 to 2018, the stock price has a clear downward trend, but the trading volume tends to be higher and there are more peaks 
during this period. This may indicate potential activities such as stock news, analyst downgrades and insider trading. 

Table 1 


Summary Statistics of the Daily Negative Log Returns from Feb. 1,1993 to Aug. 10,2018 



Mean 

Range 

Std dev 

Skewness 

Kurtosis 

Nobs 

JCP 

0.032 

(-22.582,21.962) 

2.814 

0.031 

5.561 

6430 


Table 1 summarizes the fundamental properties of daily negative log return series of JCP stock. It gives us a better view of 
the performance of JCP stock. It shows a positive average daily negative log return of 0.032, which indicates a positive 
relative loss on average. 


Table 2 

Summary Statistics of the Daily Negative Log Returns Grouped by Year 


Time 

Mean 

Range 

Std dev 

Skewness 

Kurtosis 

Nobs 

1993.1.29-1998.8.10 

-0.051 

(-5.791,8.487) 

1.588 

0.006 

1.606 

1396 

1998.8.11-2008.8.10 

0.008 

(-16.217,12.981) 

2.707 

-0.344 

3.064 

2514 

2008.8.11-2018.8.10 

0.105 

(-22.582,21.962) 

3.391 

0.179 

4.766 

2518 


Table 2 shows the summary statistics of the daily negative log returns for JCP stock that are grouped into three time periods. 
From August 11, 2008 to August 10, 2018, the average daily negative log return is 0.105 which is the highest among the 
three periods. It is apparent that the average daily negative log return increases at a high rate over time. In addition, in the 
most recent decade, the range of the daily negative log returns is the largest, from -22.582 to 21.962, with a standard 
deviation of 3.391, which indicates high volatility of the stock. 
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2.2 Test for normality 

To test for the normality of the data, we use the Q-Q (quantile-quantile) plot to see if the empirical distribution of the daily 
negative log returns is consistent with the normal distribution. 


Normal Q-Q Plot 



Theoretical Quantiles 

FIGURE 3: Q-Q plot of the daily negative log returns for JCP against normal distribution 

Fig. 3 is a Q-Q plot of the empirical distribution of the daily negative log returns (y-axis) against normal distribution (x-axis). 
It can be observed from the plot that the empirical distribution of the daily negative log returns displays heavier tails than the 
normal distribution. Therefore, normal distribution is not an ideal fit for the negative log returns. 


T Q-Q Plot (df=3) 



FIGURE 4: Q-Q plot of the daily negative log returns for JCP against Student’s t-distribution 

We also test the empirical distribution of the daily negative log returns against the Student’s t-distribution using Q-Q plot. 
From Fig. 4, we can see that with degrees of freedom of 3, the empirical distribution of the daily negative log returns has 
lighter tails than the Student’s t-distribution. It is apparent that the Student’s t-distribution is a much better fit for the negative 
log returns. 

To confirm our results, we perform Shapiro-Wilk normality test (1965) [15] to check for the normality of daily negative log 
returns for JCP. The results are shown below. 


Table 3 

Shapiro-Wilk Normality Test of Daily Negative Log Returns for JCP 


Statistic 

p-value 

Test Result 

W = 0.95186 

<2.2e-16 

The series does not come froma normal distribution 
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From Table 3, we can see that the p-value from this test is extremely small. For a significance level of 0.05, we reject the null 
hypothesis that the data comes from a normal distribution. This is consistent with our result from the Q-Q plot. 

2.3 Test for autocorrelations 

The autocorrelation coefficient (ACF) and the partial autocorrelation coefficient (PACF) are very important for us to check 
the specifications of the model used in analyzing the data. With the presence of autocorrelations, we can use autoregressive- 
moving-average (ARMA) models to fit the data. The ACF and PACF graphs for JCP negative log returns and squared log 
returns are shown in Fig. 5. From the graphs, we can see that the daily negative log return series exhibits weak 
autocorrelations, whereas the squared log return series have indication of strong autocorrelations. 
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FIGURE 5: Sample ACF and PACF for JCP daily negative log returns and squared log returns 


In order to confirm the results, we use the Ljung-Box test devised by Ljung and Box (1978) [16] to test the autocorrelations 
of the negative log return series and the squared log return series. The Ljung-Box test statistic is usually represented by 

2 

Q(m), where m is the number of lags tested. If the test statistic Q(m) is greater than X a or the p-value from the test is 

smaller than the significance level of a , we reject the null hypothesis that the data are independently distributed and 
conclude that the data exhibits autocorrelations. 

The Ljung-Box test results are displayed in Table 4. The p-values are small for daily negative log returns and extremely close 
to zero for daily squared log returns at lags of 5, 10 and 15. At significance level of 0.05, we reject the null hypothesis and 
conclude that there are strong autocorrelations within our data. 

Table 4 


Ljung-Box Test for JCP Daily Negative Log Returns and Squared Log Returns 



Daily Negative Log Returns 

Daily Squared Log Returns 

m 

X -squared 

p-value 

X -squared 

p-value 

5 

17.287 

0.003986 

279.51 

<2.2e-16 

10 

23.141 

0.01024 

395.42 

<2.2e-16 

15 

26.14 

0.03658 

501.98 

<2.2e-16 


From the above statistical tests, we can see that the negative log return series {L t } has non-normality and relatively strong 

autocorrelations. For time series with these properties, a powerful test - the ARCH test was introduced by Engle (1982) to 
evaluate the significance of ARCH effects of the data, see [4]. Here we perform the ARCH test in Table 5. We can see that 
the p-values for daily negative log return series are very small at lags of 5, 10 and 15. Therefore, we reject the null hypothesis 
at significance level of 0.05 and conclude that there are significant ARCH effects in our data. 
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Table 5 


ARCH Test for JCP Daily Negative Log Returns and Squared Negative Log Returns 



Daily Negative Log Returns 

Daily Squared Log Returns 

m 

X -squared 

p-value 

X -squared 

p-value 

5 

41.76 

<2.2e-16 

0.8852 

0.4899 

10 

25.23 

<2.2e-16 

0.5928 

0.8212 

15 

18.95 

<2.2e-16 

0.4136 

0.9761 


III. Value at Risk with GARCH and EGARCH Model 

3.1 Methodology 

ARCH model was proposed by Engle (1982) [4] to deal with the model’s time-varying volatility and heteroskedasticity of 
the errors. It was extended by Bollerslevn (1986)[6] to the generalized autoregressive conditional heteroskedasticity 
(GARCH) model. The GARCH model is used to approximate conditional variance using a linear function of the past squared 
residuals. The model still has some shortcomings because it ignores the leverage effect of return volatility. To reflect the 
asymmetry of returns, we applied the EGARCH model proposed by Nelson (1991) [7], in which the volatility can react 
asymmetrically to positive and negative returns. 

Let {L t } be the daily negative log returns of JCP stock as defined in Section 1 and {F f } be the past information about the 
return series up to time t. Since there is volatility and leptokurtosis in the data, we assume that the conditional mean of { L f } 
follows an autoregressive average model AR (1) and the conditional variance of {L t } follows an univariate GARCH model 
or EGARCH model. We represent {L t } as follow: 


(L t =/u t +cr t £ t 


the innovation {s t } are white noise process with zero mean and unit variance, and we assume of follows a GARCH-type 
model. In this paper, we assume [s t } to follow normal and Student’s t-distribution respectively. The conditional mean JU t is 
defined as: JLl t — E(L t I F t _ x ) , and the conditional variance <J t is defined as: erf = Var(L f I F t ) . 

The GARCH (p, q) model is given by: 

p q 

A = «o + H a j r llj +X^-i , forp,q >0 (3) 

7=1 i= 1 

where T] t — <J t £ t ,p is the order of GARCH and q is the order of ARCH process, a ■ and /? are parameters and we expect 
their sum to be less than 1. 

We use GARCH (1,1) model and EGARCH (1,1) model to fit the data. Petra Posedel (2005) [18], Richard A. Ashley and 
Douglas M. Patterson (2010)[19], Joel koima, Peter N Mwita, and Dankit K Nassiuma (2015) [20] have shown that the basic 
GARCH(1, 1) model is sufficiently suitable for most financial time series. From equation (3), we can easily get the equation 
for the GARCH (1,1) model: 

of =a 0 + nU + P\°)-\ ( 4 ) 

The EGARCH (p, q) model is given by the following formula: 
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p <1 

log a? = a 0 + Yj l a A-j + rj 0 Vj I -E I rj t _j I)] + £ 0 t log cr, 2 .,. , tor q > 0 (5) 

7=1 i=1 

From equation (5), we derive the EGARCH (1,1) model: 

log of = a 0 + n, , + y x (I t] l \-E\ rj t , 1) + /?, log cr , (6) 

The results of fitting data using GARCH and EGARCH models are illustrated in Table 6. The log likelihood log(L) shows 
that the serial correlations in the conditional means and variances are sufficiently explained by the specified GARCH and 

EGARCH models. The positive coefficient Y\ i n the EGARCH model implies the presence of a leverage effect. 


Table 6 

Estimation Results in GARCH and E GARCH Models for JCP 


GARCH type 

GARCH Model 

EGARCH Model 

Distribution 

Normal 

Student’s t 

Normal 

Student’s t 


-0.041362 

0.004727 

-0.000560 

0.021645 

0, 

0.012557 

0.002382 

0.007003 

-0.001397 

a o 

0.010825 

0.018960 

0.009296 

0.008382 


0.023102 

0.031519 

0.046362 

0.039226 

A 

0.975893 

0.966760 

0.996278 

0.995399 

Y\ 



0.048972 

0.065170 

l°g(L) 

-15080.79 

-14763.57 

-14995.24 

-14719.55 


3.2 Estimation of VaR 

After fitting the data with GARCH and EGARCH models, we are able to predict the volatility of the daily negative log 
returns for JCP. We then proceed to estimate the value at risk (VaR) of the stock. As introduced in Section 1, VaR is a 
popular method of measuring risk of an investment, because it is easy to interpret and clearly a relevant concept in assessing 
the risk. It represents the maximum dollar amount that is expected to be lost at a predetermined confidence level for a 

specific period of time. In this study, we calculate the VaR by the quantile method. Let {L t } be the daily negative log returns 

of JCP stock and a be the confidence level. VaR a is the a -th quantile of the distribution of the negative log returns {L f }, 
and is defined by: 

P(L t >VaR a (L t ))<l~a (7) 


Table 7 

Forecast Value at Risk in GARCH and E GARCH Models on Aug. 13th, 2018 for JCP and SPY 



GARCH type 

GARCH Model 

EGARCH Model 


Distribution 

Normal 

Student’s t 

Normal 

Student’s t 

JCP 

VaR$9 5 

4.924947 

4.355218 

5.1341 

4.523923 

VuRq99 

6.980224 

7.256845 

7.26029 

7.440851 

VaR, 

9.283975 

12.66881 

9.643528 

12.71962 

SPY 

VaR 095 

0.8465751 

0.7579894 

0.9715966 

0.8895356 

VaR, gg 

1.237037 

1.303512 

1.399687 

1.472815 

VaR, ggg 

1.674705 

2.253133 

1.879532 

2.443175 
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From the skewed property of the negative log return data, the Student’s t-distribution should fit better for the distribution of 
white noise {£ t } - To make comparisons, we use both normal and Student’s t-distributions to test the VaR. Table 7 

demonstrates the results for one-day-ahead VaR (on August 13, 2018) for JCP at quantiles 0.95, 0.99 and 0.999 respectively. 
To make comparisons, we also calculate one-day-ahead VaR (on August 13, 2018) for SPDR S&P 500 trust (SPY) [21] stock 
of which the results are shown in the same table. SPY is an exchange-traded fund (ETF) used to track the S&P 500 stock 
market index. It is the largest ETF in the world and represents all major sectors of the US market. From Table 7, it is clear 
that JCP has much higher VaR than SPY for the three quantiles. 


Series with VaR for JCP 



Series with VaR for SPY 



Figure 6: One-day-ahead VaR forecast of JCP and SPY based on the EGARCH model with Student’s t- 

distribution at quantile 95% in blue and 99% in red 


Series with VaR (95%) 



FIGURE 7: One-day-ahead VaR forecast based on the EGARCH model with Student’s t-distribution at 

quantile 95%, with JCP in red and SPY in blue 
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JCP and SPY 



FIGURE 8: Adjusted Closing Price for JCP in red and SPY in blue 


From Fig. 7, it is clear that on most of the days, JCP has much larger VaR than SPY. In addition, the VaR of JCP stock 
becomes higher and deviates further from the VaR of SPY in recent years. Fig. 8 shows a comparison graph of the adjusted 
closing price for JCP stock and SPY stock. As the price of SPY increased since 2010, the price of JCP displayed a decreasing 
trend. Combining the analysis of the VaR and the price plot of JCP, we find that JCP stock is highly risky for investment. 

IV. Calendar Effect 

A calendar effect [22] is an economic effect that seems to be associated with the calendar. If the returns of a stock vary from 
season to season, it is said to have seasonal tendencies and calendar effect. In order to investigate any calendar effect of JCP 
stock, we divide our data into four groups on a quarterly basis. Each quarter contains three months, the first quarter includes 
January to March expressed as Ql, the second quarter includes April to June expressed as Q2, and so on. 


JCP-Q1 


JCP-Q3 




JCP-Q2 


JCP-Q4 




FIGURE 9: The adjusted closing price for four quarters of JCP 


Fig.9 shows the time series plots of the daily negative log returns for the four quarters. From the graphs, we can see that the 
negative log returns of the fourth quarter Q4 are the most volatile among the four quarters. Furthermore, there are many large 
positive peaks in Q4, especially from time of 1200 and more, which indicates that in recent years, JCP stock experienced a 
lot of losses in Q4. 
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Table 9 


Summary Statistics of the Daily Negative Log Returns for JCP by Quarter 



Mean 

Range 

Std dev 

Skewness 

Kurtosis 

Observation 

Qi 

-0.068 

(-22.582,18.592) 

2.884 

-0.440 

6.851 

1579 

Q2 

-0.008 

(-16.101,21.962) 

2.689 

0.040 

6.524 

1643 

Q3 

0.124 

(-14.813,18.105) 

2.663 

0.342 

4.904 

1618 

Q4 

0.080 

(-16.217,16.674) 

3.011 

0.218 

3.888 

1590 


Table 9 shows the summary statistics of the daily negative log returns for the four quarters. From the table, we can see that in 
the third and fourth quarter Q3 and Q4, JCP stock has positive average daily negative log return with that in Q3 to be the 
largest. It indicates that JCP stock generally experiences loss in Q3 and Q4. Although the average daily negative log return in 
the first and second quarter Q1 and Q2 is negative, the magnitude is relatively small. 

After checking the normality and autocorrelations of the data for the four quarters, we fit the data of each quarter with 
GARCH and EGARCH models based on normal distribution and Student’s t-distribution respectively. Then we calculate 
one-day-ahead VaR at quantile 95% and 99% respectively of the stock negative log returns for each quarter. The results are 
shown in Table 10. 

Table 10 


Value at Risk of Four Quarters in GARCH and EGARCH Models for One-Day-Ahead Period 



GARCH type 

GARCH Model 

EGARCH Model 


Distribution 

Normal 

Student’s t 

Normal 

Student’s t 

Qi 

VaR 095 

5.5055 

5.1341 

6.0122 

5.6980 

VaR Q99 

7.8419 

8.8800 

8.5437 

9.6957 

Q2 

VaR 095 

6.3523 

5.9994 

7.2951 

6.4139 

VclRq 99 

9.0205 

10.126 

10.3406 

10.7452 

Q3 

VaR 995 

3.6436 

3.2475 

3.6474 

3.3087 

VdR 099 

5.1297 

5.5203 

5.1112 

5.5915 

Q4 

VaR Q95 

6.9963 

6.6221 

6.8536 

6.5923 

VaR 099 

9.8818 

10.9387 

9.6649 

10.7553 


From Table 10, the third quarter Q3 has the smallest VaR for both 95% and 99% quantiles, while the fourth quarter Q4 
exhibits the largest VaR for both quantiles with Student’s t assumption. The VaR for the first and second quarter Q1 and Q2 
is large as well, despite the negative average daily negative log return in Q1 and Q2 shown in Table 9. Therefore, none of the 
four quarters performs well, because none of them exhibits both positive average log return and stability. Furthermore, it is 
an interesting phenomenon that JCP has positive average relative loss and high risk in Q4, because Q4 is usually the best 
season for traditional retailers to make money. It implies that JCP is not competitive with traditional retailers, which supports 
our previous analysis that JCP has a high investment risk. 


V. Conclusion 

In this paper, we analyze the fundamental statistics and the VaR of the daily negative log returns of JCP stock. We use 
GARCH and EGARCH models to fit the data with normal distribution and Student’s t-distribution assumptions of white 
noise respectively. Based on the models, we estimate VaR of the stock at quantiles 95%, 99% and 99.9% to measure the 
riskiness of the stock. Comparing the VaR of JCP and SPY, we find that the VaR of JCP stock is much higher than SPY, and 
the price of JCP is significantly lower than SPY, especially in the past ten years. The summary statistics also indicates that 
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the daily negative log returns of JCP averaged over ten-year period increased over time and reached the highest during the 
past ten years. These results help us conclude that the risk of investing in JCP stock is very high. 

By analyzing the calendar effect of JCP, we find that although the returns in the four quarters are different, the performance 
of the stock is not good for either of the four quarters. It is surprising that Q4 has positive average daily negative log return 
and the highest one-day-ahead VaR among the four quarters. In general, for traditional retailers, the fourth quarter should be 
the most profitable, because Thanksgiving and Christmas are at the end of the year, and all goods have great discounts during 
the festivals. The fourth quarter is usually the quarter that Americans spend the most. However, even in the fourth quarter, 
J.C. Penney cannot make a profit either. It even suffered losses during Q4 and had the highest VaR in a year. It reinforces our 
previous conclusion that JCP is at immense risk for investment. 
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